AS = 2 decays of B meson in MSSM and two Higgs doublet model 
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In view of the extreme smallness of AS — 2 transitions of B meson in the Standard Model, we consider their 
occurence in several extensions of it. Thus, we analyze the three - body — + K^K^tt'^ and two - body 
B- K*-K*°, B- K-K°, B- K*-K°, ^ R-R*" decay modes both in the Standard Model and in 
the Minimal Supersymmetric Model with and without TZ parity conservation and in two Higgs doublet models. 
All five modes are found to have a branching ratio of the order of 10~^^ in the Standard Model, while the expected 
branching ratio in the different extensions vary between 10~^ — 10~^, for a given reasonable choice of parameters. 



The rare B meson decays are very important 
in current searches of physics beyond Standard 
Model (SM) |l[]. Recently, it has been suggested 
|2|-Q to investigate effects of new physics possi- 
bly arising from b ssd or 6 ^ dds decays. 
As shown in Ref. the b — > ssd transition 

is mediated in the standard model by the box- 
diagram and its calculation results in a branching 
ratio of nearly 10""'^^, the exact value depending 
on the relative unknown phase between t, c con- 
tributions in the box. The authors of Refs. 
have calculated the b ssd transition in various 
extensions of the SM. It appears that for certain 
plausible values of the parameters, this decay may 
proceed with a branching ratio of 10~^ — 10"^ 
in the minimal supersymmetric standard model 
(MSSM) and in two Higgs doublet models [|). 

Thus, decays related to the b — > ssd tran- 
sition which was calculated to be very rare in 
the Standard Model, provide a good opportu- 
nity for investigating beyond the Standard Model 
physics. Moreover, when one considers super- 
symmetric models with 7?,-parity violating cou- 
plings, it turned out that the existing bounds 
on the involved couplings of the superpotential 
did not provide any constraint on the b — > ssd 
mode Recently, the OPAL collaboration 
has set lower bounds on these couplings from the 
establishment of an upper limit for the — ^ 
K-K-TT+ decay BR{B~ K-K-tt+) < 8.8 x 
10^^. Here we investigate the branching ratio 



of this decay mode in MSSM, with and without 
TZ parity and two Higgs doublet models as pos- 
sible alternatives to the SM. Then we comment 
on another possibilty for the observation of the 
b — > ssd transition: the two body decays of B^ . 
First, we proceed to describe the framework used 
in our analysis in which we concentrate on MSSM, 
with and without TZ parity and two Higgs doublet 
models. 

The minimal supersymmetric extension of the 
Standard Model leads to the following effective 
Hamiltonian describing the b ssd transition 
ii 

n = CMSSM{s"f''dL)isjf,bL). (1) 

where we have denoted 

Cmssm = - "^Fl" m^hix) + 66/6(x)] (2) 
216™ 7 
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with X = rn^/m?-, and the functions f^ix) 

and fQ{x) are given in The couplings df^ 

parametrize the mixing between the down-type 
left-handed squarks. At the scale of 6 quark mass 
and by taking the existing upper limits on 5^ 
from ^ and |^ the coupling Cmss m is estimated 
to be \Cmssm\ < 1-2 X 10^^ GeV^'^ for an aver- 
age squark mass toj = 500 GeV and x = 8, which 
leads to an inclusive branching ratio for b ssd 
of 2 X 10^^ The corresponding factor cal- 
culated in SM, taking numerical values from |^] 



and neglecting the CKM phases is estimated to 
be ICsmI ^4x10-12 §. 

The authors of |^ have also investigated be- 
yond MSSM cases by including R- parity violat- 
ing interactions. The part of the superpotential 
which is relevant here is W = X[jf.LiQjdk, where 
i,j,k are indices for the families and L,Q,d are 
superfields for the lepton doublet, the quark dou- 
blet, and the down- type quark singlet, respec- 
tively. Following notations of ||] and the tree 
level effective Hamiltonian is 

^ = - E ^K32K2li-^BbL){sLdR) 

n ^ri 

+ K2lK32i-SRdL){sLbR)]. (3) 

The QCD corrections were found to be important 
for this transition For our purpose it suffices 
to follow 1^ retaining the leading order QCD re- 
sult fqcD ^ 2, for mo = 100 GeV. 

Most recently an upper bound on the spe- 
cific combination of couplings entering (|^) has 
been obtained by OPAL from a search for the 
B~ ^ K-K-7r+ decay § EJIA^^gaA^Tai P + 
|A'„2iA'„*39p)^^^ < 10"''. Here we take the order 



n32l 

of magnitude, while the OPAL result is 5.9 x 10~^ 
based on a rough estimate r{B^ K~ K~it~^) ~ 

1/4 r{b ssd). 

The decay b — > ssd has been investigated us- 
ing two Higgs doublet models (THDM) as well 
These authors found that the charged Higgs 
box contribution in MSSM is negligible. On the 
other hand, THDM involving several neutral Hig- 
gses |l^ could have a more sizable contribution 
to these modes. The part of the effective Hamil- 
tonian relevant in our case is the tree diagram 
exchanging the neutral Higgs bosons h (scalar) 
and A (pseudoscalar) (see for details P,p^). 

1 1 _ _ 
'Hth = T;isb^sdi^{sd){sb) 

2 mf^ 



-^(S75C?)(S75&)), 



(4) 



The coupling defined in Q as a Yukawa cou- 
pling of the FCNC transitions di ^ dj . In our es- 
timation we use the bound \£,sb£,sd\/'m'jf > 10~^° 
GeV^^, H = h,A, which was obtained in |^ by 
using the AriiK limit on S.bd/'m-H and assuming 



lU/mnl > 10-3. 

We proceed now to study the effect of Hamilto- 
nians (||), (||) on the various two body AS" = 2 de- 
cays of charged B - mesons. In order to calculate 
the matrix elements of the operators appearing 
in the effective Hamiltonian, we use the factor- 
ization approximation Jll|-[l3t , which requires the 
knowledge of the matrix elements of the current 
operators or the density operators. Here we use 
the standard form factor representation [|l2[|ll|l of 



the matrix elements described in detail in y,|li 

For the Fi and Fq form factors appearing in 
the decomposition of the matrix element of the 
weak current between two pseudoscalar states, 
one usually assumes pole dominance |l^,|l^. For 
the vector and axial vector form factor, appear- 
ing in the decomposition of the matrix element of 
the week current between the vector and pseu- 
doscalar states, we use again pole dominance 
[12,14|. The relevant parametrs are taken from 
|ll|,|l|l F(f^'(0) = 0.38, A^^*{0) = 0.32. For 
the calculations of the density operators we use 
derivatives of the vector or axial-vector currents 
|l5|| . In our numerical calculation we use /k = 
0.162 GeV, gK' = 0.196 GeV^ |l|. Now we turn 
to the analysis of the specific modes. 

We denote O = (s7^(l - 75)^) (57^(1 - 
75)5), and then we use H = CO with C being 
l/'iCMSSM, l/4CsM [|l5|. Using factorization 
and introducing s = (ps — /ci)^, t — (ps — ^2)^ = 
and u = {pb — P-kY , one finds for the —>■ 
K^K^TT^ decay 

{K-{k^)K-{k2)^+{j>^)\0\B-{pB)) = 
Fi" {s)F^^ {s)[ml + 2m^ - s - 2<- 



)] 



+Fo^-(,s)F-«(.) 



BK, 



(5) 



Within MSSM the branching ratio is found to be 



BR{B- 



MSSM 



< 4.7 X 10" 



(6) 



while SM gives this rate to be 5.2 x lO^^''. The 
MSSM which includes TZ parity breaking terms 
can occur in this decay. The matrix element of 



the operator O-jz 
found to be 



{K-{h)K~ik2)7r+{p^)\On\B-{pB)) 



(7) 



(ms - md){mi, - nis) 

Taking the values of the quark masses as in 
rUb = 4.88 GeF, = 122 MeV, rud = 7.6 MeV 
and using the bound given above, we estimate 
the upper hmit of the branching ratio BR{B~ 
K-K-Ti+)n < 1.8 X 10"^. This hmit can be 
raised to 6 x 10~^ for the upper bound on the 
coupHngs of 5.9 x lO"'^ given in Q. 

The two Higgs doublet model, with the limit 
l^sb^sdl/mjj > 10" ^° GeV~'^, results in a branch- 
ing ratio of the order 10" 

The long distance effects (LD) are usually sup- 
pressed in the B meson decays. However, in any 
search of new physics one has to include their 
contributions also [Q. In the case oi B~ — > 
K~K~Tr^ decay, we have analyzed two contri- 
butions 01 : (I) the box diagram, which is essen- 
tially the LD analog of the SD calculation in the 
standard model |^ of the b — > ssd transition. 
(II) the contribution of virtual and "ttO" 

mesons, via the chain B" ^ ii'""£>0"("7rO") 
— > K~ K~7r'^ . This contribution arises as a se- 
quence of two AS* = 1 transitions and may lead to 
final K~K~7T~^ state as well. It is therefore nec- 
essary to have an estimate of its relevance vis - a - 
vis the "direct" AS* = 2 transition. The box dia- 
grams contributes to the real and imaginary part 
of the amplitude for the B~ — > K~K~tt^ decay. 
We have found that the real part, for a reason- 
able value of the cut-off parameter A ~ 10 GeV 
results in the rate 8 x 10"^^, while the imaginary 
part of the amplitude leads to the rate 6 x 10"^^. 
The largest nonresonant contribution of the " 



or "TT"" exchange comes from the "£)°", givmg a 
branching ratio smaller than 10"^^. Therefore, 
we have shown that the long - distance contribu- 
tions to B~ — > K~K~'K^ in the SM are smaller or 
comparable to the short - distance box diagram, 
and have the branching ratio in a 10"^^ — 10"^^ 
range. This is a most welcome feature since it 
strengthens the suitabilty of the S" K^K^tt^ 
decay as an ideal testing ground for physics be- 



yond the standard model, as originally suggested 
in ref. [|. 

We briefly discuss the two - body AS* = 2 de- 
cays of B" meson. Although in principle two 
body decays would appear to be simpler to an- 
alyze, there is the complication of K'^ — mix- 
ing. Hence one needs also a good estimate for the 
b — > ssd transitions as well. For the analysis of 
pseudoscalar meson decay to two vector mesons 
B" K*^K*^ it is convenient to use helicity for- 
malism (see details in ||l|]). Within MSSM model 
the branching ratio becomes < 6.2 x 10"^, while 
SM gives this rate to be 6.8 x lO-^"*. The 7^ - par- 
ity term described by the effective Hamiltonian 
(4) cannot be seen in this decay mode when fac- 
torization approach is used, since the density op- 
erator matrix element {K*^\{sd)\Q) vanishes. The 
two Higgs doublet model also cannot be tested in 
this mode due to the same reason. 
The use of factorization technique described 
above gives following results in the case of B~ 
K*^K° decay: within MSSM the branching ra- 
tio is straightforwardly found to be BR{B~ 
K*-K°)mssm < 1-6 X 10"9 (l|], which is com- 
parable to the SM prediction of Ref. ||l|] for 
the A5 B" -> K*-K° decay given as 
BR{B- K*-K°) = lxl0"^ 5xl0"^ 2x10"^ 
obtained for the number of colours Nc — 2, 
Nc = 3, Nc — oo, respectively. The SM calcula- 
tion for the AS — 2 transition leads to BR{B^ 
K*-K°)sM = 1.7 x 10-^^. The MSSM which in- 
cludes TZ parity breaking terms can occur in this 
decay. The estimation of the upper limit of the 
branching ratio gives BR{B^ K*^ K")-jz < 
4.4 x 10"*. This limit can be raised to 1.5 x 10"^ 
for the upper bound on the couplings of 5.9 x 10"^ 
given m |]. The two Higgs doublet model (H) 
gives for the limit \(,sb^sd\/m'jj > 10"^" GeV^, 
a branching ratio of the order 10"^^. Due to spe- 
cific combination of the products of the scalar 
(pseudoscalar) densities this is the only decay 
which has nonvanishing amplitude within the fac- 
torization assumption. 

For the B~ K~K*'^ decay mode the branch- 
ing ratio in MSSM is constrained to be BR{B~ — > 
K*'^)mssm < 5.9 X 10"^ in comparison with 
SM result 6.5 x 10"^**. The ampHtude calculated 



in MSSM including 7^ breaking and THDM van- 
ishes, due to vanishing of the matrix element of 
the density operator for K*^ state. 
The B~ K~K^ decay offers the following: 
the branching ratio for MSSM is found to be 
BR[B'^ K-K°)mssm < 2.3 x 10"^, in com- 
parison with the 2.5 x 10"" found in the SM. The 
matrix element of the R parity breaking MSSM 
operator has nonvanishing value and the con- 
straint on the coupling constants < 10^'* gives 
the bound 9.4 x 10"^, while for the bound of 
5.9 X 10^"' for the coupling constants the rate 
BR{B- ^ K-K°)n can reach 3.3 x IQ-^. 
One might wonder if the long distance effects 
are important in two - body AS" = 2 B~ de- 
cays. We have estimated the tree level contribu- 
tion of the D{D*) which then goes into K{K*) 
via weak annihilation. We found that these con- 
tributions give a branching ratio of the order 
10"^* and therefore they can be safely neglected. 
One might think that the exchange of two in- 
termediate states D{D*), K{K*) can introduce 
certain long distance contributions. In decay 
B "D" "isT" ^ "if" "if" the first weak vertex 
arises from the decay i? — » " D" " K" and the sec- 
ond weak vertex (see e.g. |^) can be generally ob- 
tained from the three body decays of D ^ KKK . 
Therefore, we are quite confident to suggest that 
the long distance effects are not important in the 
two - body AS' = 2 i? decays. 

We can summarize that in the B^ — > 
K-K-TT+ and two - body B' decays, the MSSM 
with the chosen set of parameters gives rates of 
the order 10~^ — 10~*, while the TZ parity breaking 
terms in the MSSM can be seen only in the B~ — » 
K-K-TT+, B- K*-K° and K' de- 

cays. Let us turn now to the possibility of detect- 
ing these decay modes. The B~ K~K~Tr~^ 
seems to be the best candidate, since the other 
modes we discussed, have a K'^ in the final states 
which complicates the possibilty of a detection 
because of —K^ mixing These are the 

modes which as we mentioned are more difficult 
on the experimental side. The THDM model 
can give nonvanishing contribution in the case 
of B- K-K-TT+ and B- K*- K° decays, 
with a rate too small to be seen. Thus, we con- 



clude that the B K K 7r+ decay is an ideal 
candidate to look for physics beyond SM and that 
the B- K*-K*°, B-^ K-K*° decays offer 
this possibility also. 
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